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Abstract. A generalisation of the Faulhaber polynomials and Bernoulli num- 
bers related to elliptic curves is introduced and investigated. This is applied 
to compute the density of states for the classical Lame operators. 

1. Introduction 

Johann Faulhaber, " Rechenmeister der Stadt Ulm", in his book |H] pubhshed in 
1631 announced a remarkable fact that the sums of the odd powers of the first n 
natural numbers can be expressed as the polynomials of their sum. Namely, if we 
denote 5*1 = 1 + 2 + ... + n (n^ + n)/2 as A then 

53 = l^ + 23 + ...+n3 = A2, 

5*5 = 1^ + 2^ + ... + = (4A=' - A2)/3, 

^7 = 1^ + 2^ + ... + = (12A* - SA^ + 2X^)/6 

and in general 

02m-l = 1 +2 +...+n = -t,n[X) 

for some polynomials Fm , which are called Faulhaber polynomials: 
i^i(A) = A, F2(A) = A2, F3(A) = iA2(4A-l), F^iX) ^ ^X\6X^ - 4X + 1), 

F5(X) = iA2(16A3-20A2 + 12A-3), Fe(X) = iA2(16A'*-32A2 + 34A2-20A + 5). 
5 3 

In full generality this claim was first proved by Jacobi in 1834 jj^. For a nice 

discussion of this story and the effective ways to compute the Faulhaber polynomials 

we refer to a very interesting paper jl7) by Donald Knuth. 

For the general k the sums of powers Sk-i{n) = 1*^^^ + 2*^^^ + • • • + n'^~^ can 

be expressed through the classical Bernoulli polynomials (see e.g. 

= {Bk{n + l)-Bk)/k, 

where Bk = Bk (0) are the Bernoulli numbers. Bernoulli polynomials can be defined 
through the generating function: 

ze"'-" _ Bk{x) 
- 1 ^ k\ 

k=0 

and have the symmetry property Bk{l — x) — {—l)''Bk{x). Faulhaber's claim is a 
simple corollary of this symmetry. In fact, Faulhaber polynomials are related to 
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the even Bernoulli polynomials in a simple way: 

B2mix + 1) = 2mF^ (^T^) + ^2m- 

Faulhaber did not discover Bernoulli numbers but his work was known to Jacob 
Bernoulli and cited in his famous treatise in Ars Conjectandi in 1713. 

Recently, D. Fairlie and one of the authors |S] discovered an interesting relation 
of the Faulhaber and Bernoulli polynomials with the theory of the Korteweg-de 
Vries (KdV) equation 

Ut - 6uUx + Uxxx = 0. 

It is known since 1967 due to Gardner, Kruskal, Miura and Zabusky [El EDI that 
this equation has infinitely many conservation laws of the form 

(1) Im[u]^ J Tm{u,Ux,Uxx, ■■■,Um-2)dx, 

where Tm[u\ are some polynomials of the function u and its ^-derivatives up to 
order m — 2 : 

Ti=u, ^2 = M^ ^3 = ^il^+2u^ T4 = U2'^ + 10uui^ + 5u'^... 

They are uniquely defined by some homogeneity property modulo adding a total 
derivative and multiplication by a constant (see the next section). This constant 
can be fixed by demanding that Tfn{u, u^, u^x, Um-2) = pl^^ ^ function of 

derivatives of order less than m — 2. 

The KdV equation has famous solutions known as solitons, the simplest of which 
is a one-soliton solution 

u = -2sech^(x - At), 

corresponding to the initial profile u{x, 0) = — 2sech^a;. The main result of ^ is the 
following formula, relating the Faulhaber polynomials Fm with the integrals of the 
KdV equation: 

22m 

/™[-2Asech\] = (-1)"- -FmiX). 

2m — 1 

In the present paper we introduce a new class of polynomials, which we call 
elliptic Faulhaber polynomials. They are related to the periodic traveling waves 
(known also as cnoidal waves) of the KdV equation: 

u{x, t) = 2p*(x — ct), 

where p*(x) satisfies the differential equation 

(2) (p'J^ = 4p^ -5iP* -52P* -53 

with gi = —c. The function p, (x) differs from the classical Weierstrass elliptic 
p-function by adding a constant. 

The elliptic Faulhaber polynomials are defined as 

(3) .F™(A|r,7) = /T„[2Ap*(z)]dz, 

where the integral here is taken over a cycle 7 on the corresponding elliptic curve 
F given by the equation — AX^ — giX"^ — — with X — p^{z),Y — p'^{z). 
Note that since the integrand has all the residues zero, one can consider 7 here 
as an element of the first homology group i?i(r, Z) and the integral is a linear 
function on it. Equivalently, one can say that the integrand is a differential of a 
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second kind and thus, by the de Rham theorem, determines a special element of 
the first cohomology group H^{T, C), polynomially dependent on the parameter A. 
Let 

be the corresponding periods of the basic second-kind differentials on this curve. 
The variables gi, (72, <73, w, ^ can be considered as the coordinates on the space 
U ~ {(r,7)}, so Tm = J'mi^', 91, 92, 93,1^,0- Define the weight of the variables 
A, 91,92, 93, t^, £. as 0, 2, 4, 6, — f , -l-f respectively. 

Theorem 1. The elliptic Faulhaher polynomial J>„ is polynomial in all variables 
\,gi,g2,93,uj,^ with rational coefficients, homogeneous with weight 2m + 1. When 
92 = ga = it reduces (up to a factor) to the classical Faulhaher polynomial 

(4) (A; 31 , 0, 0, c^, = - T^^9T''^F„^ (A). 

2m — f 

Here are the first 4 elliptic Faulhaher polynomials: 

^3 - (-^.9?C + ^.9i52cc)A2(4A - f ) - ^g2^X\3X - 2) + ^g3u;X'{2X - 3), 
•^4 - i^gU - ^9h2io) X\&X^ - 4A + f ) - ^5152? A2(26A2 - 29A + 9) 

+ -gi93UjX^{iX^ - 2A - 3) + —glujX^{2hX^ - 40A + 24)] - —g^£^X^{bX^ - 15A + 9). 

If we consider the Weierstrass reduction gi — Q we come to the reduced version 
of these polynomials {X; g2, g3,uj,v) = •?^m(A; 0, 32, ffs, 77), where p* = p, ^ 
reduces to = p{z)dz = ({ui; g2, ga), p and ( are the standard Weierstrass 

elliptic functions [22! ■ One can see the explicit form of the first 8 of them in the 
Appendix. 

These polynomials have double zero at zero (for to > 1) with the second deriva- 
tive at zero proportional to 

which we call elliptic Bernoulli numbers. Note that in contrast to the so-called 
BernouUi-Hurwitz numbers BH2k |12[I15| . which are functions of only 32,53, these 
numbers B2m — B2m(,g2, 53, w, 77) depend also on uj and rj: 

K K 3 2 R 2 2 36 

5 5 7 7 

On the discriminant (when two of the roots of the equation AX^ — g2X — 53 = 
collide) they reduce up to a simple factor to the usual Bernoulli numbers: 

>32n{^9l, ^5?, ^,^+^91^) = -B2ngit 

Our main motivation for considering these polynomials came from the theory of 
the classical Lame equation 

-y" + n{n + l)p{z)y = Ey, 
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well known since XIX century (see 0] EHl)- The classical procedure going back 
to Hermite and Halphen allows one to find the general solution y{z) for given 
n = 1,2,3.... (see [21 El ^oi the recent results in this direction). Our approach is 
somewhat different: we want to say something about the Lame equation for all 
values of the parameter n. 

More precisely, we consider the density of states p{E) for the corresponding Lame 
operator L = —-j^+n{n+l)p{x). The density of states is one of the most important 
spectral characteristics of the Schrodinger operator with periodic potential (see e.g. 
|23j). From the finite-gap theory [131 121| it is known that the density of states of 
the Lame operator with integer n has the form 

for certain polynomials P„, i?2„+i. On the other hand, the high-energy asymptotics 
for p can be given in terms of the KdV integrals (see \21\). which in this case are 
our (reduced) elliptic Faulhaber polynomials: 

piE)='[l + ^t'-^^]. A^!^. 
2ttVe 22fc-i E'' ^' 2 

Let p — —-^ be the average of the Weierstrass p- function over real period T = 2uj. 
We have the following 

Theorem 2. The coefficient = ak{n) in the numerator Pn{E) = E"" -\- aiE""^^ + 
a2i?"^^ + • • • + a,i of the Lame density of states |31) is a polynomial in n of de- 
gree [^]. Its coefficients are polynomials with rational coefficients of g2,g3 and p, 
homogeneous of weight 2fc, where the weights of g2, gs are 2, 4, 6 respectively. 

The procedure is actually effective and gives the explicit form of the first coeffi- 
cients afc(n) for all n: 

nin + 1) _ 

02 = -^(n- l)n(n+l)(6 + 25n+16n^), 

03 = (n-2)(n-l)n(n+l)(--S7T(45-h243n+247n2+64n3)-^n(n-M)(27+16n)) 

3360 960 

More of these formulas, which we believe to be new, are given in section 6 of this 
paper. 

In the last section we discuss also some open problems and conjectures about 
the elliptic Faulhaber polynomials. 

2. The densities of KdV integrals 

The KdV integrals Ik play a fundamental role in the theory of the one-dimensional 
Schrodinger operator L = —-j^ + u{x). They appear in the asymptotic expansion 
of the trace of the resolvent of i? = {pi — L)^^ (see Gelfand-Dikii (T]) or, in the 
periodic case, of the density of states or quasi-momentum (see Novikov et al fH\ 
and section 6 below) . In this section we will discuss the algebraic properties of the 
corresponding densities Tfc mainly following the classical paper by Miura, Gardner 
and Kruskal (2111 ■ We should warn the reader that our u differs from the variable 
in that paper by a constant factor, so the arithmetic of the coefficients is slightly 
different. 
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The densities of the KdV integrals Tk — Tk{uo,ui, ...,Uk~2) are polynomials in 
u = uq and its derivatives Ux — ui,Uxx — U2,--- up to order k — 2. In general 
they are defined up to a total derivative, which does not affect the integrals. To 
eliminate this freedom, let us introduce the notion of canonical form ,20;. Consider 
the following grading: for the monomial Mq"°u!|"^ ...u™'' we define its rank as 

j=o 

In other words, the rank of u is 2 and each derivative gives additional +1. 

A form T{uq, ui, Ur) is called canonical if it is rank-homogeneous and is irre- 
ducible in the sense that there are no monomials u^qu]^...u]^ in T with the highest 
derivative appearing linearly (i.e. ik = !)• 

It is clear that if the highest derivative Uk appears linearly one can always elim- 
inate such a term by doing integration by parts (which is equivalent to adding a 
total derivative). 

Theorem 3. (Miura, Gardner and Kruskal 2Q:J For any k = 1,2,.... the KdV 
equation has a non-trivial integral Ik with a density Tk — Tk{uo,ui, ...,Uk-2), 
which is a rank-homogeneous polynomial of rank 2k. The canonical form of the 
corresponding density Tk is defined uniquely up to a constant factor. 

The highest derivative term Uk~2 of the corresponding Tk is it^ ^- We normalise 
these densities by choosing the coefficient at this term to be 1. It turns out that 
this makes all the other coefficients integer, which will be important for us. 

Theorem 4. All the coefficients of the KdV density Tk in the normalised canonical 
form are integers. 

Although this fact was probably known to the experts in KdV theory, it did not 
play any role before. We have not found the proof in the literature although it is 
not difficult. 

Proof of theorem 4- One can define the KdV densities by the following recurrent 
formula: 

k—m—l 

(6) (Tm+i = -0-„ — ^ CTkCTyn-k 

fe=l 

with ai = u (see e.g. PTl). These densities are not in the canonical form but are 
rank-homogeneous and obviously have integer coefficients. The even cr are known 
to be total derivatives [2], so only odd a lead to non-trivial integrals. The claim 
is that 

(7) Tk = {-If-^o^k-i, 

where the polynomial cr2fc-i irreducible equivalent of 02k-i- We should only 

show that the coefficient at u\_2 in Tk is equal to I. 

From we have (T2fe-i — ~'^2k-2 ~ Y^'^i=t <^i'^2k-2~i- The total deriva- 
tive term cr2j._2 vanishes when we take the irreducible form o'2k-i does 
not contribute to the term u^ g- The only contributions to it^_2 come therefore 
from — X]i=i'^ criO'2fe-2-i- Since the only term of degree I in ai is (— I)*~^Mi_i 
the contribution to u'l_2 comes from — ^*~J'^^'^(— l)*^^Ui_i(— I)^'^^'^^*W2fe-3-i — 
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~'"ft-2 ~ 2X]i=o UiU2k-i-i- After integration by parts UiU2k-i-i is reduced to 
(— l)'^~*~^w|_2, SO the quadratic term u|_2 in <T2k-i coefficient — 1 — 

2(-l)''E'=S~^(-l)' = (-1)''"^ foi' all k. The theorem is proved. 
Below are the first KdV densities (in the canonical form): 

Ti=u 

Ts = ui^ + 2u^ 
Ta = U2^ + lOuui^ + bu^ 
U = us'^ + 14mm2^ + TOu^Mi^ + 14?/ 
Tg = U4^ - 20w2^ + 18mm3^ - 35iti'^ + 126u^U2^ + ^20u^ui^ + 42tt^ 

The following proposition gives an explicit formula for the last coefficient in these 
polynomials. 

Proposition 1. The coefficient oJuq in the canonical form is ^1(^22)1 for k> 1 
and 1 for k = l. 

Proof. The term Uq appears with exactly the same coefficient (up to a sign) in 
T/c as in a2k-i since the integration by parts does not affect it. Consider again 
the relation a2k-i = ~'^'2k-2 ~ (^i(^2k-2-i- The total derivative term (j'2k-2 
does not contain any term Uq. The only contributions to Uq come therefore from 
— X^-Zi*^^ (7i(72k-2-i- This leads to the following recurrence relation for the coef- 
ficients bk in Tfc = bkUQ + : 

i=k-l 

bk= bibk-i 

i=l 

with bi = 1. For the generating function G{t) = J2k=T ^kt'^ this means that G{t)'^ — 
G{t) + t = Q, i.e. G{t) = ^~^2^~^* . Expanding G{t) as a power series and equating 
the coefficients, we obtain the required expression for bk- 

3. The elliptic Faulhaber polynomials 

Consider now the special potential of the form m = 2Ap* , where — p*{z;gi,g2,g3) 
is the elliptic function satisfying the differential equation 

{p'*f = 4p* - gipl - g2p* - 53- 

This function differs from the classical Weierstrass function p{z;g2-,gi) with the 
equation 

{p'Y = ^P^ - 92p- h, 
where 52 = 52 + and 53 = 53 + i^5i52 + ^9i, by adding a constant: 

P*{z;9i,g2,g3) = p{z;g2.h) + 

However it will be convenient for us to keep the additional parameter gi and consider 
the elliptic curve F in the non-reduced form 

y2 = 4X3 _ ^^^2 _ g^-^ _ 



ELLIPTIC FAULHABER POLYNOMIALS AND LAME DENSITIES OF STATES 



7 



Choose any cycle 7 on F, which does not pass through the poles of the function 
< and define the elliptic Faulhaber polynomials as the integrals 



(8) ^rn(A|r,7) = (j, T„,[2Xp*{z)]dz. 

Since the integrand has all the residues zero, this integral can be considered as a 
linear function on the first homology group of F or equivalcntly as an element of 
the first cohomology group H^{r, C). Thus, the elliptic Faulhaber polynomials can 
be considered as A-dependent sections of the canonical cohomology bundle over the 
space of the elliptic curves £ = {F}. 

We choose the following basis in the first cohomology of F : 

(9) ^^nfdz = -(b—, £_ = -- (b p {z)dz = -- 



2 J 2/y'" 2 J ° ' ' 2 J Y 

Now we are going to prove our Theorem 1 claiming that J>„ ~ T„i{\, ffi, (?2i 53, ^, 
are polynomial in all variables with rational coefficients. 

Proof of Theorem 1. We will need the following lemma about the derivatives 
of the function p,, which can be easily proved by induction using the differential 
equation Q for the function p» and its corollary 

(10) p" = &pI - giP* - ^32- 

Lemma 1. The derivatives of the elliptic function p* have the form 

(11) pl^''^ = yl^(p,;gi,52,53), pl^''^^' = ^^'(p,; 51, 32, 53)p* 

for some polynomials A^, with rational coefficients. The polynomials A^, satisfy the 
following recurrence relation 

(12) ^fc+i = (4p^ - gipl - 52P* - g3)Al" + {6pl - 5ip, - ^52)^^', 
where the derivative of A^, is taken with respect to p* . 



We know that the KdV density Tm{uo, Um-2) is a sum of monomials 
with rank 

j=fe 

r = y~^(2 + i)mj = 2m. 
3=0 

From lemma 1 it follows that if we substitute in such a monomial u = 2\p^^,{x) we 
will have an expression of the form Q(p*, 32, <?3)A^^, where M = X]j=o "^i ^^"^ 
Q is a polynomial in p*,5i,<72,33 with rational coefficients. Indeed, the fact that 
the rank of Tk is even means that the first derivative of p* appears only in even 
powers and thus can be eliminated using the equation It is easy to see also 
that if we define the grading of p*, (71,32,53 to be 2,2,4 and 6 respectively then Q 
is homogeneous of degree 2m. 

Let us consider now the integrals 



mu 



= f p"iz)dz. 

Proposition 2. The integrals K* satisfy the following recurrence relation of the 
third order 

(8n - A)K: = (2n - 2)g,K:_, + (2n - 3)g2K_, + {2n - 4)g3K:_, 



8 



M-P. GROSSET AND A. P. VESELOV 



with initial terms Kq — 2uj, K* = —2^, K2 = ^92^ — ■|5iC- 

Indeed, using ^ we have K* = / p'^-^{z)pl{z)dz = § p"~2(z)(^p"(z) + f p* + 
%)dz - ^/p»-i(z)d^ + a^§p-^'\z)dz + \§pT\z)p'i{z)dz. Integrating by 
parts the third integral gives § p^~'^{z)p'l{z)dz = —{n — '2)§ p"^'^(z)(p'^(z))^dz — 
-{n-2)§pr\z) (4p3 (z) - 51 pI [z) - 52 p* [z) - g3)dz ^ -{n-2) [AK* - g,K*_, - 
g2K*-2 - g3K*-3)- Therefore K* = fK*_, + I^K*_, - ^{AK* - g^K*_, - 

g2K:^2 - 93K:_,) = -^K* + V;5i/C-i + ^92K*-2 + ^53i^;-3, which 
leads to the required recurrence relation. The expression for Kq,Ki are obvious; 
the form of A'2 follows from the relation (|10|l . 

Corollary 1. The integrals K* have the form 

(13) = a1")(5i, 52,53V -M"^(.9i, 52, 53)C, 

where yli™^(5i, 52, 53), -B*'"'' (51, 52, 53) are polynomials with positive rational coeffi- 
cients. 



Here are the next few integrals: 

_ I , „ , 1 

^3 - 30 ^^-^-^ ■ 30 



K; = ^(5152 + 653)^ - ^(25? + 952)e, 
1 o o , 1 



Kl - ^m92 + 25gi + mgm)^ - ^(35^ + 265152 + 6053)^, 

^5 = ^(4g?g2+335i52+245?53+1685253)^-^(85?+1025252+14752'+3005i53)^. 

Combining all this we have the proof of the first part of the theorem. To prove 
the second part let us consider the case 52 — 53 = 0. The equation ^ then becomes 
simply 

(p'j2 = 4p3-5ip2 

with the solution p*(z) = jj^j^i^^, where a is related to 51 by 51 = —Aa . 

For the vanishing cycle the integral £, — —■^ § p^:{z)dz is identically zero while 
the integral £, — —-^ p*{z)dz taken over the real line and corresponding to the 
diverging period lo has the finite value a. In fact, it is convenient to shift the real line 
in the complex plane by i^, then the function becomes p*(x + — — a^sech^aa;, 
which (up to a coefficient ^) is the soliton profile (see the Introduction). 

Now we can apply the result by Fairlie and one of the authors 5 , which claims 
that 

o27n 

(14) /™[-2Asech2x] = (-1)"- -F^{\). 

2m — 1 

Let us remind here the idea of the proof. Consider the solution of the KdV equation 

Ut - GuUj: + Uxxx = 

with initial data 0) = —n{n + l)sech^a;. It is well-known that for integer n 



i{x, t) 2 ^ fc^sech^fc(a: - Ak^t - Xk) 



as < ^ 00 for some constants Xk (see e.g. 0, page 79). Now we can use the 
fact that the integrals are actually the conserved quantities, so /m[u(a;,i)] — 
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Im[u{x, 0)]. Because these integrals have the property Im[a'^u[ax)\ — a^™ ^Im[u{x)\ 
this immediately gives us the equaUty 

n 

Im[-n{n + l)sech2a;] = /,„ [-2sech\] ^ /fc^m-i^ 

k=l 

Now to derive the formula (|14|) we need only to show that/™ [-2sech^a;] = (_i)™22", 
1), which can be done in various ways. 

All this leads to the formula l@J, which completes the proof of Theorem 1. In 
fact, one can prove a slightly stronger result: 

Proposition 3. The elliptic Faulhaber polynomial has the form 

where F„i{X) is the usual Faulhaber polynomial and the dots mean the terms of 
lower order in gi. 

There are two more interesting special cases: lemniscatic when ffi = 53 = and 
equianharmonic when gi — g2 — 0. Both of them are specialisations of the reduced 
version of the elliptic Bernoulli polynomials, which we are going to look at in more 
detail in the next section. 

4. The reduced elliptic Faulhaber polynomials. 
The reduced form of the elliptic Faulhaber polynomials corresponds to gi = : 

In that case the function becomes simply the standard Weierstrass elliptic func- 
tion p(z] g2, 93), so by definition the reduced elliptic Faulhaber polynomials are 

(15) F^i\\T,j)= <fT^[2\piz)]dz. 

where 7 is any cycle on the elliptic curve F given by the algebraic relation 

y2 = 4X3-52X-.g3, 

which does not pass through the poles of the p-function. 

Remark. One can consider also another reduction replacing the Weierstrass 
function p by the Jacobi elliptic function fc^sn^, appearing in the Jacobi form of 
the Lame equation |29| . The corresponding Jacobi reduced version of the elliptic 
Bernoulli polynomials can be defined by setting g^ — : 

Fi{\;9i,92,i^,Cj = •?^m(A;5i,.g2,0,cj,^). 

It corresponds to the elliptic curves F with a chosen half-period and has the ad- 
vantage that one can easily see the hyperbolic limit (52 — 0). However we prefer 
the Weierstrass version as it is more canonical and standard in the mathematical 
literature. 

Theorem 5. The elliptic Faulhaber polynomials F^ have the form 

= A„^{kg2,g3)i^ + S,„(A; 52, 53)'?, 
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where the sum is taken over all non negative integers k, I satisfying 2k + SI = m 
and 2k + 31 = m — 1 respectively and j4{^™-'(A), b'"^^\X) are some polynomials in A 
of degree m with rational coefficients and having double zero at \ = Q. 

The explicit form of the first 8 polynomials can be found in the Appendix. We 
are going to look in more detail at the highest and the lowest degree coefficients in 
A of the polynomials . 

The highest coefficient of is proportional to Kn — § p{z)^dz. Setting gi = 
in the proposition 2 from the previous section we have 

Proposition 4. The integrals Kn satisfy the following 2 term recurrence relation 
of third order 

{8n - 4)X„ - {2n - 3).g2i^„-2 + (2n - 4)g3i^„-3 
with the initial data Ko = 2ljj, Ki = —2ri, K2 = ^g2^- They have the form 

(16) /C™ (32, 53^-5^™^ (52, 53)^, 

where A^''"'\g2, gs), B^"^\g2, gz) are polynomials with positive rational coefficients. 

There is no explicit way to solve this recurrence (except for the lemniscatic and 
equianharmonic cases, see below), so we would like to discuss here another way 
to compute Kn, going back to the classical work by Halphen JI] (see chapter VII 
page 203). 

It is based on the following relation between powers of p and its even derivatives: 



r=0 ^ ' 

where B^"^ are coefficients given by the Halphen recurrence relation (see I2(jp 
(18) 

(„) _ (27i-2r-2)(2n-2r-l) 2n - 3 („-2) , "-2 («-3) 
(2n-2)(2n-l) +4(2n-l)^'-2 2(2n - 1) 33 

with n > 0, r = 0, ...,n and with Br"-* ~ for r < or r > n, Bq"-* = 1 and 
Sj") = , for any n. By construction B^"^ are polynomials in g2 , gs with rational 
positive coefficients. We will call them Halphen coefficients; see the first of them in 
the Table 1. 



Table 1 . The Halphen coefficients 



n 





1 


2 


3 


4 


5 


6 


r=0 


1 


1 


1 


1 


1 


1 


1 


r=l 























r=2 








1^52 








^52 


r=3 











To93 




¥% 


^53 


r=4 














336^2 


240^2 


gf)0 52 


r==5 

















Mff2.93 


15^^'^? 2 


r=6 




















492852 + 55^3 
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This leads to the following expression for Kn in terms of Halphen coefficients: 
(19) if„ = / p"{z)dz = 2B(")a; - 25^1^7,. 



Indeed, since the complete integrals of the total derivatives are zero, the only non- 
zero contribution comes from the constant and p terms in the right-hand side of 
the relation (|17|l . This is of course related to the general de Rham theorem, which 
implies that on the elliptic curve T any differential p"(z)(iz is cohomological to a 
linear combination of dz and p{z)dz. 

(n) (n) 

Note that these particular Halphen coefficients -B„ and satisfy the same 

recurrence relation Q as Kn'. 

2n - 3 n~2 

but with initial conditions q_i = 0, ao = 1,q;i = for a„ = si""* and a_i = 
0, Qfo = 0, ai = I for a„ — iJ^lV 

In the lemniscatic case when 53 = this recurrence equation can be easily solved 
and leads to the following formula for Kn- 

(20) if„ = 2^ l[{ik-2r+lfg[^\{l~r)u;^2rrj). 

Similarly in the equianharmonic case when 172 — 0, we have three cases depending 
on the congruence of n modulo 3: 



n — 3fc + I [|.] 



2lf]-i 11 2n-6fc-F5 

k—1 



93 



/orn^l[3], K^^-^fl 



n-Sk + 1 i^] 

93 V, 



2n - 6fc + 5 

k—1 

forn=2 [3], ^ 0. 

In the next section we analyse the coefficients of the elliptic Faulhaber polyno- 
mials at the lowest degree A^, which can be considered as an elliptic generalisation 
of the Bernoulli numbers. 



5. The elliptic Bernoulli numbers 

We define the elliptic Bernoulli numbers B2n = ^82)1(52, 53, ry) as the following 
integrals 

(21) B2n = ^j{j—^p{z] 92.93)? dz. 

They are simply related to the lowest degree coefficient at in the reduced elliptic 
Faulhaber polynomial F^(A), which is equal to 8K2m-2- Note that the coefficient 
at lowest degree A^ in J>n(A) for m > 1 can be obtained by replacing the coefficients 
52, 53, -n in S2m-2 by 52 + ^ff? and 53 + ]^gi52 + ^9i, S. + respectively. 
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The terminology is justified by the fact that on the discriminant, where two roots 
of the polynomial AX^ — g2X — gs collide, these numbers reduce to the classical 
(even) Bernoulli numbers: 

R - ^ B - ^ R - ^ R - ^ R - ^ R - ^ R - 
-fc>0 — i, i>2 — 7, r>4 — — — , i36 — -77;^ — ^TTTT) ^10 — 7:77, ^12 — 



6' 30' 42' " 30' 66' 2730' 

Recall that all the odd Bernoulli numbers except Bi = —1/2 are zero. 

Proposition 5. For n > 1 the specialisation of the elliptic Bernoulli numbers on 
the discriminant is 

(22) B2n{j^gl, ^5?, u;, ^ + ^g,u;) = -B^n 5?^, 
where B2n cltc the usual Bernoulli numbers. 

This follows from Theorem 1 and the relation between the coefficients of 
and J-jn mentioned above. One can show this also more directly using the integral 
formula for the Bernoulli numbers found in [5] . 

Let us look now at the properties of the elliptic Bernoulli numbers in the general 
case. 

Proposition 6. The elliptic Bernoulli numbers satisfy the following recurrence 
relation: 

= (-l)"-^(2n-l)![<+^)c.-(i3(/'+i)-i3("))r;]-(2n-l)!^ ^ ^^ ""[^ g2n-2., 

r=2 " ^ 

where S^"'' are the Halphen coefficients. 

Indeed, integrating by parts, we have §{p{z)^"~^^)'^dx = (-1)""^ / p{z)'''^"''^^ p{z)dz, 

so 

(23) B2n = f p{z)(^--^^p{z)dz. 
Re-arranging the formula p7|) gives 

= (2n- l)![p" - ^^|£^p(--6) _ ... _ :^p(2) _ B^%p-Blr^ 

Multiplying the above formula by p{z) and taking the contour integral over a 
cycle gives 
(24) 

B2n^ {-lT-'^^^^^^[j p-+\z)dz-j B^-^p{z)dz]-{2n-\)^^^^ 

Using the relation K^+i = § p'^+^{z)dz = 25^"^^^ - '^B^n'^'^^ri and § p{z)dz = 
—277, obtain the proposition. 

The recurrence relation in this proposition can be considered as an elliptic ana- 
logue of the recurrence relation for the usual Bernoulli numbers: 

3=0 V / 

where ('"^^) are the binomial coefficients. 
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There is another way to compute the eUiptic BernouUi numbers if we assume 
that we know aheady the integrals Kn- Let us use the following very convenient 
notation: for a polynomial D{x) = X^jLo '^j^"' ^^^"^ ^ sequence Kq, Ki, . . . define 
D[K] as follows 

N 

D[K]=Y,d,K,. 

Let Dn — x^„_i(a;; 33) = 0, 52, 53), where A^. are the polynomials 

defined in Lcmma^by the recurrence relation (|12|l . Then from (|23ll we immediately 
have 

Proposition 7. The elliptic Bernoulli numbers can be expressed through Kn as 
follows: 

Here are the first few elliptic Bernoulli numbers: 

„ 3 2 „ 2 2 36 

^2=3^2^^, B4 = --gsuj + -g2r], Bq = -g2UJ - —gsT], 
5 5 7 7 

36 24 2 „ 72, 3 , „^ 2160 
Bs = — -929?.^ + -^92V, Bw = —{92 + 1833)^ ^5253??, 

298512 2592 

612 ^3233^+^(49(7f+750g|)r,, 614 = 21652(gi+365|)t^-9072g253^, 

S16 = -i^ff3(1039g| + 450032)^^ + 1^32(5395! + 180005f)r7. 

The special lemniscatic (173 = 0) and equianharmonic {92 = 0) cases can be seen 
straight from here, but we can not present a closed formula for Bn for general n. 

Remark. In the classical case the coefficients of the Faulhaber polynomi- 
als Fm{X) can be computed from the fact that in the corresponding polynomial 

Fm{^) = 2k(^2m(x + 1) " ^2™) = X^™"! + ^{B2m{x) - B2„), all the CO- 

efficients at odd powers of x except x^"^~^ are zero (see ^T]). This is related to 
the fact that the coefficients of the Bernoulli polynomials are proportional to the 
Bernoulli numbers and all the odd Bernoulli numbers except Bi = —1/2 are zero. 
A simple check shows that this is not true in the elliptic case: already for JFg the 
coefficient at x^ in jTg (2-^) is ~ 2^30; (which is of course zero if 92 = 93 = 
but not in general). This means that the elliptic Bernoulli numbers do not play the 
same role in the elliptic case as they do in the usual one. 

It is worthy to mention that there is another elliptic generalisation of Bernoulli 
numbers: the so-called Bernoulli- Hurwitz numbers BH2k |12[ I15L I22 | . They are 
defined as 

Bi/2fc -2fc(2fc-2)!cfc, 

where Cfc are the coefficients of the Laurent series of the Weierstrass p-function at 
zero : 

P(^) = ^ + C2Z^ + C3Z^ + ... 

These coefficients satisfy the recurrence relation 

3 ''"^ 
= (2fc+l)(A:-3) ^ 
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for fc > 4 with C2 = 92/20,03 = 33/28 (see EH)- 

The BernouUi-Hurwitz numbers are related to Halphen coefRcients in the follow- 
ing way. Recall that the principal part of a Laurent series is the sum of the terms 
with negative powers. Since the principle part of p*^^™) is ; we have from 

the relation H17|l the following 

Proposition 8. The principal part of the Laurent series p" — {-^ + C2z'^ + c^z'^ + 
)" is 

n-l (n) 
Df 



y2n 



r-=0 

Since (J^ + C2z2+C3z4 + )" = ^ + ^ + ^ + ^^ - -^ 

ncrj+n(n-i)c2C3 ^ ^^^^^ allows US to cxprcss the BernouUi-Hurwitz numbers through 

Halphen coefficients and vice versa: 

B^") = nc2 = ng2l2Q ^^BH^, 

S(")=nc3- 7133/28 =^i3i/6, 

r(")_„, , ^('^-1) „2_ n{3n-l) _ n{3n-l) 

t>. — ncA H Co — Ca — orifi, 

* * 2 ^ 2 11520 

B^"^ = nc5 + n{n - l)c2C3 = -(22n - 19)c5 = "(^^"-19) ^^ 
5 o -r V y 2 3 g V ; o 1209600 

Note that these expressions for B). are only valid when r < n. 

6. Application: the density of states of the Lame operators 

The Lame operator L„ = — ^ + n{?2 + l)p{z) has been a classical object of 
investigation since the XlX-th century although its spectral properties were appre- 
ciated much later. In 1940 Ince [TJ proved that the Lame operator (on the real 
line shifted by an imaginary half-period) for integer n has exactly n gaps in its 
spectrum. The ends of the gaps are given by the zeroes of certain polynomials 
R2n+i{E) , which we will call spectral. The problem of finding these polynomials 
was investigated in the classical works by Hermite and Halphen (see I^Hl) but still 
continues to be of interest (see e.g. recent papers by Belokolos and Enolski [2] and 
Takemura |^|2S] and references therein). 

In this section we will discuss the density of states for the Lame operators, which 
is one of the most important notions in the spectral theory of the Schrodinger 
operator with periodic potential (see e.g. |23j). In the finite-gap case, in particular 
for the Lame operator i„, the density of states has the form 

(25) piE)^l- , 

where Pn{E) = E^ + aiE"~^ + a2i?"^^ + ... + a„ for some coefficients ai, 02, a„ 
and R2n+i{E) = ]XZo"{E — Ei) = J2k=o"^^ 6fc_E^"+^~'^ is the spectral polynomial. 
We are going to explain how to find the coefRcients Uk as a function of n in terms of 
the elliptic Faulhaber polynomials and the coefRcients of the spectral polynomials, 
which were investigated in jlOI. In particular we prove the Theorem 2 from the 
Introduction. 
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Let US start with the known resuhs about the density of states for the Schrodinger 

,2 

operator L = + u{x) on R with periodic potential u{x) of period T (see e.g. 

|23[ 12 Ip. The spectrum of the operator L is known to be continuous and has a 
band structure with in general an infinite number of gaps. Consider the restriction 
L^-^') of L on {—R,R), for R > 0, with Dirichlet boundary conditions at ±R. The 
spectrum of L^^' is then discrete. Denote the eigenvalues as £„(L^^)) and define 
the integrated density of states N(E) as the limit 

where tt{£'„(L*^-'^)} is the counting function of the discrete spectrum of L^-^-\ It 
is known that such limit exists and has the following asymptotic expansion as 

where Ik[u\ — JrpTk{u,uu2, ■ ■ ■ ,Uk-2)dx are the KdV integrals evaluated over 
period T. 

The density of states is defined as the derivative of the integrated density of 
states: 

It is known (see 23 ) that for the periodic potentials u{x) the integrated density of 
states is related in a very simple way to the so-called quasi-momentum p{E): 

(27) N{E) = 

TT 

Recall that the quasi-momentum p{E) appears naturally in relation with the Bloch 
(Floquet) solutions of the Schrodinger equation Lip = Eip, which have the form 
4'Eix) = exp'^^-^)"^ 4>ix)i 4>ix) being a periodic function of period T. 

If u{x) is a finite gap potential with n gaps then the derivative of the quasi- 
momentum p{E) has the form j21|: 

dpjE) ^ 1 PnjE) 
dE 2^R^-^) 

where Pn{E) — E" + aiE"^^ + a2-E"~^ -I- ... -l-a„ for some coefficients oi, a2, a„ 
and R2n+i{E) = ni=o(-^ ~ ^0 where Eq, Ei, E2n are the end points of the 
spectrum intervals of L. 

Thus in that case, the density of states has the following form 

(28) p{E) = ^—M^)=. 

^R2n+l{E) 

Note that from (|26|) it follows that it has the following expansion at infinity 

(29) PiE)--^ + -^ypk-l)''^''^ 



2TTy/E ttT/B^ '{4E)''' 

Consider the case of the Lame operator L„ = + n(n + l)p. The integrals 

Ik[u] are then the elliptic Faulhaber polynomials F^(X), where A = g^j^j 
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T = 2{jj. Thus, the high-energy asymptotics for the densities of states of the Lame 
operator can be given in terms of the elhptic Faulhaber polynomials as follows: 



1 1 °° 



2^/B 2tj^^22fc-i 

Thus in principle, if we know all the Faulhaber polynomials then we know the 
density of states of Lame operators L„ /or all n. In particular, one can use a 
general procedure based on the theory of continued fractions going back to Stieltjes 
to reconstruct it. However it does not look effective enough, so we will take a 
different approach. 

First, we express the coefficients hu of the spectral polynomials using the so- 
called elliptic Bernoulli polynomials ^U]. More precisely, in JO] we have proved 
that the coefficients bk of the Lame spectral polynomial 

2n 

R2n+l{E) = 1[{E - E,) = F^n+l ^ ^^£;2„ ^ ^^g^n-l ^ ^ ^^^^^ 

are polynomials in n, 172, 93 with rational coefficients. The degree of b^ in n is [^] : 

61 = 



&2 = -]^"(" + 1)(2" - l)(2n + l)(2n -f 3) 

63 = -Z77;n{n + l)(2n - 3)(2n - l)(2n + l)(2n + 3)(2?i + 5) 
840 

64 = ^2 „(n-l)(n+l)(2n-l)(2n+l)(2n + 3)(56n^+76n^-94n2 + 201n + 630). 
201600 

Re-writing the high energy asymptotics of the density of states of the Lame 
operator i„ in the form 

1 Pn{E) _ 1 ri , "1 , / 1 . ai&2 63 X 1 , 

and equating the coefficients term by term with the expansion ( 131(1 gives 

F^{\) 7]^ r]n{n + l) 
ai = — = A = , 

Aw Ul Ll) Z 



04 



03 



b2 


3 


F2^(A) 


0,2 = 

2 


+ 16 




0162 


|. 


5 F3^(A) 


~ 2 




26 uj ' 


362 


0163 


64 7 [A] 




2 


2 28 oj 



2 

From this we find recursively 

n{n + 1) 
"1 



02 = — —{n - l)n{n + 1)(6 + 25n + 16^2), 
480 
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aa = --||-(n-2)(n-l)ri(n+l)(45+243n+247n2+64n=^)-^(n-2)(n-l)n2(n+l)2(27+16n), 
ooDU you 

04 = „„5^,^^ (n-3)(n-2)(n-l)n(n+l)(-2520-12942n-10315nV4565n^+6880n^+1792n^) 



3225600' 



3)(n - 2)(n-l)n2(n+ 1)2(600 + 563n+ 128^2), 



13440 

where p = — 5 average of p(a;) over a period. Note that the coefRcient is 

divisible by {n + l)n{n — l)(rt — 2)...{n — fc + 1), which is related to the fact that 
L„ has n gaps for integer n. All this implies Theorem 2. 

Remark. The differential p(E)dE — tt — f"^'^-''^'^ has a simple algebro-geometric 

meaning: it is an Abelian differential with a second order pole at infinity normalised 
by the condition that all the periods over gaps (6-periods) are zero (see e.g. Novikov 
et al PU). This means that once the equation of the spectral curve is known one 
can find the coefficients of the numerator by solving a linear system of equations 
with the coefficients being some standard hyperelliptic integrals on this curve. From 
our results it follows that these coefficients actually can be expressed polynomially 
in terms of the standard elliptic integrals and the parameters of the initial elliptic 
curve. This is of course related to the reduction problem for the Lame spectral 
curves, see e.g. pi^. 

One can see the explicit form of the first 7 coefficients and the first 5 Lame 
densities of states in the Appendix. They are in a good agreement with Belokolos- 
Enolski and Takemura calculations |21 who used a modified version of the 
classical approach to this problem. 

7. Some open questions and conjectures. 

The coefficients of the classical Faulhaber polynomials are known to have alter- 
nating signs (see ^31). We conjecture that the same is true for the reduced elliptic 
Faulhaber polynomials. More precisely, recall that in the Weierstrass form these 
polynomials have the form = Am{X; (?2, ffs)^ + -Bm(A; 52, ffs)??, where 

Conjecture 1. The coefficients of the polynomials A^j^^^ {X) and B^j^\X) have al- 
ternating signs. 

For the first 8 polynomials one can check this from the explicit form given in the 
Appendix. We have checked that a similar property holds for the Jacobi reduced 
versions (but not for the general case, see e.g. Ta in the Introduction). An ideal 
proof would be to find a combinatorial interpretation for these coefficients (cf. 
Gessel and Viennot 8 and Knuth J3 i^i the classical case). It might be easier to 
prove it first for the elliptic Bernoulli numbers (cf. Proposition 4 above). 

Conjecture 2. The elliptic Bernoulli numbers B2m have the form 

B2.n = (-l)"-i(i('"n52,53)«^ - B'-"^\92.9^)ri). 
where the polynomials A(™^((72, 53), B^"^\g2, 53) have positive rational coefficients. 

Another interesting question is the behaviour of the elliptic Faulhaber polyno- 
mials on a real line. Actually it is more instructive to look at the corresponding 
polynomials J-mi ^ t^ ) as functions of x rather than A — ^ . Indeed, we know 
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Figure 1. Left: ^six) for the usual normalised Faulhaber poly- 
nomial, Right: $g^(a;) in the elliptic lemniscatic case 



that the integer values x — n should play a special role here, being related to the 
" finite-gap" values of the parameter in the Lame equation. Recall also that in the 
usual case these polynomials coincide with ■^{B2m{x -f 1) — B2m), so their graphs 
are simply shifted graphs of the corresponding Bernoulli polynomials. 

Conjecture 3. As m tends to oo for x in a finite interval on the real line 

l__cos27ra 

In the hyperbolic limit 32 = 53 = (i.e. for the usual Bernoulli polynomials) 
it is known to be true (see e.g. EHl)- The normalisation constant 2K2m-2 is 
chosen to guarantee the correct second derivative (chosen to be 1) at zero. Another 
justification for this conjecture is that the asymptotic expansion ( I31|l should be 
convergent for integer n (due to the finite-gapness property), which means that the 
numbers J^m{ " 2 " ) should be relatively very small. 

In the Fig. 1 we show the graphs of the normalised polynomials 



«'m(a;) = 



2S2m-2 

for m — 8 in the usual and lemniscatic elliptic case gi — gs = 0. 

The elliptic Faulhaber polynomials have further generalisation related to the 
so-called Treibich- Verdier potentials 

3 

u = mo(mo + l)p{x) + ^ mi{mi + l)p{x - uji), 
1=1 

where uji are the half-periods of the Weierstrass elliptic p-function. The corre- 
sponding KdV integrals Ik[u] are polynomials in four variables — mi(mi + l), i — 
0, 1, 2, 3 with the coefficients depending on ei, 62, 63, w, rj. We are planning to discuss 
them elsewhere. 

We believe also that the elliptic Faulhaber polynomials and elliptic Bernoulli 
numbers are interesting from a number-theoretic point of view and would like to 
mention the papers ^1 El |221 in this relation. 
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9. Appendices 

9.1. Appendix A. The first 8 reduced elliptic Faulhaber polynomials 

The lemniscatic and equianharmonic cases can be read off by putting 53 = and 
52 = respectively. 



pw ^ - r] [AX] 



- g^n [fA2(3A-2)] 
^ ~ + 9^^ [|A2(2A-3)] 



+ gluj [^A2(25A2-40A + 24)] 
^ " - [f A2(5A2-15A + 9)] 



- ah 



[^A2(49A3 - 140A2 + 168A - 72)] 
+ [i|A2(28A3 - lOSA^ + 126A - 54)] 



+ 92^ [^^^(45A* - 200A=^ + 416A2 - 400A + 144)] 
pw ^ - g2g3V [Wa2(87A* - 515A3 + llTgA^ - 1206A + 450)] 

+ glco [||A2(56A*-420A=^ + 1197A2- 1368A + 540)] 



- 92V [i|A2(847A5-5390A^ + 17248A3-30536A2 + 26928A-9072)] 

pW ^ + g^gau [3^A2(9526A5 - 71995A^ + 250404A3 - 472428A2 + 433224A - 149256)] 

- 93V [f^A2(220A5-2310A^ + 10395A3- 22770A2 + 22572A-8100)] 

+ giuj [|i A^(1521A® - 13104A^ + 59696X^ - 165568A^ + 269568A^ - 224640A + 72576)] 
pW ^ - glgzr} [||A2(2171A^ - 22477A^ + 113295A'' - 336492A3 + 570492A2 - 485784A + 158760)] 

8 

+ g^gloJ [f A^(182A® - 2184A^ + 12285A* - 38844A^ + 68094A^ - 58968A + 19440)] 

9.2. Appendix B. The list of the first 7 coefficients ak{n) of the numerator 
Pn{E) of the density of the Lame states 

Here we use the notations P = —-^^ Uk{n) = (n + l)n(n — 1) . . .{n — k + 
1), ak{n) = ak{n)/Uk{n). 

ai = \p, 

a2 = -^(6 + 25n+16n2), 
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as = - gflo (45 + 243n + 24771^ + 64n^) - f|§n(n + 1)(27 + 16n), 

a4 = 32^00 (-2520 - 12942n - lOSlSn^ + A565n^ + 6880n'' + 1792n^) - -^^n{n + 
l)(600 + 563n+128n2), 

as = i77|g|no (-28350 - 145305n - 98919n2 + 130400n3 + 185250n'' + 7848071^ + 



^2 - 



11264n6) + 645^1200 "(" + l)(-22050 - 19707n + 321771^ + 7328n3 + 1792n^), 

= 3228825600 (585728n^ + 6077568n/' + 240557107?;'^ + 42381080^4 + 22989372n-'' - 

21506058rt^-26135595n-4677750)- gg4^^^^3^ooo (4100096n^+23905024n^+14017296n^- 
192219241n5-520562096n4-391295859n3+180468864n2+310981356n+62868960)+ 
Y§§i^n{n + 1) (22528n5 + 171950n4 + 42704571^ + 1\hl\hn^ - 568458n - 612360), 

«7 = - i5498f 6288000 (164 3 8471^ + 160552704n^ + 520180864n^ + 153752039n*^ - 

2909673459n-^ - 6672918014^4 - 4259587899n3 + 2522284551n2 + 3574728990n + 
681080400) + g4^^4^3Q„ n(n+l)(585728n6+6709056nH29129901n4+54097586T?.-'^+ 

19647288n2-62676225n-60031125)-^3283^^4ooon(n+l)(41000967i^+25442560nS+ 
8937936n5-279814135n4-748894946n3-399303225n2+770746914n+806818320). 

9.3. Appendix C. The explicit form of the first 5 Lame densities of states 

• For n = 1 the potential is u{x) = 2p{x). 

The density of states is well-known in this case: 

n=l 1 E+p 1 E+r 



PLame{E) = 



27r ^E-i -fE+f 27r ^(E + ei){E + e2){E + es) 
For n = 2 the potential is u{x) = 6p{x), the density of states is 
„=2 E^ + 3pE-lg, 



T,^AE^ - 21g2E^ - 2753^2 + 27giE + 8I5253 
For n = 3 the potential is u{x) = 12p(a;) and the density of states is 
E' + 6pE^-fg2E-ifgs-fg2P 



2n^E^ - fg^E^ - ^^g,E^ + ^-^glE^ + '-^g^g^E^ - ^(fff - 27gl)E 

• For n = 4 the potential is u{x) = 20p(.T). The density of states is 

„=4 ,r^,_ E' + WpE^ ^g2E^ C-^gs + ^g2P)E + ^.gj - 875g3p 

PLame\^) — ; 

27T^J{E3 - 52g2E - 56O53) UlzliE^ - ^OekE - 35el - 7^2) 

• For n = 5, the potential is u{x) = 30p(a;) and the density of states is 

p5 -1- 1 c;,;:;p4 631 _ p3 / 6615 „. , 4815 „. ,-M7'2 i ^ 18117 „2 42525 „. ,-^7 i 178605 i 13365 „2,: 

-fc + ^g2-fc - (—^93 + —^g2p)li' +(—§—92 —g3p)ii^-\ § — 9293-1 —92S' 

2n,J (£2 _ 27ff|) UlzliE^ + 15efc£;2 + (3i5e2 _ 13252)-^ - 675el - 54O53) 
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